The Small World of Seismic Events 
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The understanding of long-distance relations between seismic activities has for long been of interest 
to seismologists and geologists. In this paper we have used data from the world-wide earthquake 
catalog for the period between 1972 and 2011, to generate a network of sites around the world 
for earthquakes with magnitude m > 4.5 in the Richter scale. After the network construction, we 
have analyzed the results under two viewpoints. Firstly, in contrast to previous works, which have 
considered just small areas, we showed that the best fitting for networks of seismic events is not 
a pure power law, but a power law with exponential cutoff. We also have found that the global 
network presents small-world properties. Secondly, we have found that the time intervals between 
successive earthquakes have a cumulative probability distribution explained by the non-extensive 
statistical mechanics. The implications of our results are discussed. 



INTRODUCTION 

The general belief in seismic theory is that the relation- 
ship between events that are located far apart is hard to 
be understood. However we live today in a world where 
data is being collected on most aspects of our lives and 
better yet, computer power is cheaply available for an- 
alyzing the data. When we put together the computer 
power and the data we open a series of possibilities to 
look for patterns in the data. The work on seismic data 
analysis is no different; we have now large collections of 
millions of seismic events from around the world which 
deserves analysis. In this paper we demonstrate that we 
live in a small world when it comes to these relations 
between seismic events. An event in a particular geo- 
graphical site appears to be related to many other sites 
around the world and not only other events at nearby 
sites. 

The ability to find useful information from data is not 
new and is commonly known as Data Mining. However, 
since the work from Barabasi and Albert [T] researchers 
have turned their attention not on mining the data itself 
but rather organizing the data in a network which cap- 
tures relationships between pieces of data and only then 
mining the network structure and hence the relations be- 
tween pieces of data. The network can review informa- 
tion that could not possibly be seen from mining the 
raw pieces of data. The use of networks as a framework 
for the understanding of natural phenomena is nowadays 
called Network Sciences or Complex Networks. 

In the last few years, some analysis related seismic phe- 
nomena demonstrate that earthquakes display features 



explained from the perspective of non-extensive statis- 
tical mechanics [2H5]- These networks display complex 
features that can be better understood statistically using 
the Tsallis entropy [6 . Through the analysis the spatial 
distances and the time intervals between successive earth- 
quakes using non-extensive statistical mechanics, the au- 
thors have found that two successive earthquakes are in- 
divisibly correlated, no matter how spatially far away 
they are each other [7|. 

In line with the successive earthquake model men- 
tioned above, recent studies '8, Oj have applied concepts 
of complex networks to study the relationship between 
seismic events. In these studies, networks of geographi- 
cal sites are constructed by choosing a region of the world 
(e.g. Iran, California) and its respective earthquake cat- 
alog. The region is then divided into small cubic cells, 
where a cell will become a node of the network if an 
earthquake occurred therein. Two different cells will be 
connected by a directed edge when two successive earth- 
quakes occurred in these respective cells. If two earth- 
quakes occur in the same cell we have a loop, i.e., the 
cell is connected to itself. Fig. [l] depicts a toy example 
of a network being formed. This method of describing 
the complexity of seismic phenomena has found that, at 
least for some regions, the common features of complex 
networks (e.g. scale-free, small-world) are present. How- 
ever, in spite of the importance of the results which show 
that seismic networks for some specific regions present 
small-world effects, these results are somewhat expected, 
since makes sense that areas located geographically near 
each other, are correlated. 

In this paper we have used data from the world-wide 



earthquake catalog for the period between 1972 and 2011, 
to generate a network of sites around the world. Since 
only seismic events with m > 4.5 are recorded for any 
location in the world, we then consider these significant 
events and used them in our analysis. The results were 
analyzed under two viewpoints. The first, under the per- 
spective of the complex networks theories, and the second 
using non-extensive statistical mechanics. 



THEORETICAL BACKGROUND 

Complex Netw^orks Features 

Scale-free networks are defined as one in which the de- 
gree distribution of the nodes (or vertices) in the network 
follows a power law, that is, the probability that a net- 
work will have nodes of degree k, denoted by -P(fc) is 
given by 



P(fc) 



(1) 



where 7 is a positive constant. 

Eq.[l] states that scale-free networks have a very small 
number of highly-connected nodes (called hubs) and a 
large number of nodes with low connectivity. These net- 
works exist in contrast with general random networks in 
which where the probability distribution follows a Pois- 
son distribution for networks with the number of nodes 
very large as in 
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where N is the number of nodes in the network and each 
node has an average of {k) connections, p represents 
the probability of an edge to be present in the network 
and can be shown to be approximately {k)/N. Random 
networks have nice properties but the truth of the matter 
is that most real networks are not random. 

The definition of a small-world network is yet to be for- 
malized. One of the best approaches for defining small- 
world networks is based on the work of Watts [TU] which 
states that in small-world networks, every node is "close" 
to every other node in the network. It is generally agreed 
that "close" refers to the average path length in the net- 
work, i, which has the same order of magnitude as the 
logarithm of the number of nodes, i.e.. 
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In addition, and what makes small world even more 
interesting, is the fact that these networks have a high 
degree of clustering representing a transitivity in the rela- 
tion of nodes; if a node i has two connections, the theory 
argues that the two connections are also likely to know 



each other. More formally, the clustering coefficient, Ci, 
of that node is given by: 
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where, A(i) is the number of the directed triangles 
formed by i with its neighbors and Aau{i) is the number 
of all possible triangles that i could form with its neigh- 
bors; the clustering coefficient of the entire network, C, 
is just the average of all Ci over the number of nodes 
in the network, N. In random networks the clustering 
coefficient can be estimated using the closed form 
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where (k) is the average degree in the random network. 
Last, one needs to understand the relation of these 
two characteristic to the world of seismic events. If a 
network of seismic events contains hubs, one can argue 
that the distribution of earthquakes per country around 
should also follow a power law. On the other hand, if 
the network of seismic events has small-world properties 
one can argue that there is some indication of long-range 
relations between far-apart earthquake sites. 



Summary on Non-extensive Statistical Mechanics 

Nonextensive statistical mechanics is a theory that 
explains many physical systems where the traditional 
Boltzmann-Gibbs statistical mechanics is insufficient. 
The nonextensive statistics can explain a variety of 
complex systems with characteristics such as long-range 
interaction between its elements, long-range temporal 
memory, fractal evolution of phase space, and certain 
kinds of energy dissipation. These systems have a nonex- 
tensive entropy, i.e., the entropy of the whole system can 
be greater or smaller than the sum of the entropies of its 
parts. In these cases, we use Tsallis entropy [S], which 
is a generalization of the Boltzmann-Gibbs entropy (de- 
fined later in this paper). The Tsallis entropy is defined 
as: 



S, = K 
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and W is the total number of possible configurations, pi is 
the associated probabilities, _ftr is a conventional positive 
constant, and q is the entropic index that can be under- 
stood as a measure of the nonextensivity of the system. 



This entropy violates the additivity property, i.e. if we 
have two systems statistically independent A and B, 



S,iA + B) = S,{A) + S,{B) 



K 



S,iA)S,{B), (8) 



where we can see that the quantity |1 — f?] gives the mea- 
sure of nonextensivity of the system under investigation. 
Taking the limit g — > 1 in Eqs. [6] and [8J we get re- 
spectively the standard Boltzamann-Gibbs entropy, S — 
—K^■^-^^Pi\npi, and the additivity (extensivity) prop- 
erty. It is also known that applying the maximum en- 
tropy principle to the Tsallis entropy, the probability dis- 
tribution obtained has a g-exponential form |llj . where 
the g-exponential function, eq{x), is defined by, 
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(9) 
The inverse function of the g-exponential is the q- 
logarithmic function. 
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where, for the limit g — > 1 in Eqs. |9] and [TO] we get 
back the standard exponential and logarithmic functions, 
respectively. 



A GEOGRAPHICAL NETWORK FROM SEISMIC 

EVENTS 

The use of networks to understand phenomena associ- 
ated with geographical locations has been used in many 
instances in science including diseases [H], scientific col- 
laborations [13 m], and organ transplantation [TS] to 
mention but a few. Seismic activity intrinsically linked 
to geography because todays instruments can pinpoint 
with great accuracy the location in the globe where each 
seismic event take place. 

It is important to precisely locate the geographical lo- 
cation of a seismic event but if we want to understand 
relations between events we should concentrate on creat- 
ing a network in which locations are liked based on an 
acceptable criteria. In this paper we use the same pro- 
cedure employed by [5] in their studies of earthquakes 
in specific regions of the world. The construction of the 
network is as follows. We first have to decide on what 
should represent by the nodes. Obviously our first choice 
are the sites where the earthquake took place. The prob- 
lem of doing this is that an earthquake epicenter is rarely 
located exactly in the same location and given the accu- 
racy of today's instruments we would have an infinitely 
large number of possible sites. We decide instead to de- 
fine nodes representing a larger region of the world we 
here call a cell. A cell will become a node of the network 
if an earthquake has its epicenter therein. The creation of 



edges follows a temporal order of seismic events. For in- 
stance, if an earthquake occurs in a cell Ci and the next 
earthquake in a cell C2, we assume a relation between 
Ci and C2 and we represent by the event by a directed 
edge in the network. The process continues linking cells 
according to the temporal order. Fig. [T] depicts the pro- 
cess used to create the network from seismic events. It 
is worth noting that if two successive earthquakes occurs 
in the same cell, this node will be connected to itself via 
a self-edge or a loop. 

Time-ordered sequence of seismic events 




FIG. 1. This figure shows how the network of seismic events 
is created. At the top of the figure we see a sequence of time- 
ordered seismic events. Since each event has an epicenter E 
with location (Oe, 
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to calculate 



bs) we can use Equation 
which cell in the map to be used as a node in the network. 
The nodes are linked based on the sequence of events shown 
at the top of the picture. 



The degree of the nodes (the total number of its con- 
nections) is not affected by the direction of the network. 
The nature of the way the network is constructed means 
that for each node in the network, its in-degree is equal to 
its out-degree (the exceptions are only the first and last 
sites in the sequence of seismic events but for all practical 
purposes we can disregard this small difference). 

Although the use of temporal ordering of events is not 
new in our paper, there are two main differences between 
our study and others. Firstly and most importantly, the 
region considered in our investigation is the entire globe, 
instead just some specific geographical subareas of the 
globe; this is the first worldwide study of seismic events 
using networks and consequently the first one to inves- 
tigate the possibility of long-range links between seismic 
events. Secondly, we have used a two dimensional model 
which the depth dimension of the earthquake epicenter 
is not considered, since we are interested in looking for 



spatial connections between different regions around the 
world and besides 82% of the earthquakes, in our dataset, 
have their hypocenters in a depth less or equal to 100 km. 

Before we divide the globe into cells, we need to choose 
the size of such cells particularly because we are dealing 
with the entire globe; if the cells are too small we will 
not have any useful information in the network, if the 
cells are large we lose information due to the grouping 
of events into a single network node. There are no rules 
to define the sizes. Therefore we have taken three differ- 
ent sizes, the same sizes used in previous studies [8l I16j , 
where the authors conducted studies about earthquake 
networks using data from California, Chile and Japan. 
The quadratic cells have, 5 km x 5 km, 10 km x 10 km 
and 20 km x 20 km. To set up cells around the globe, we 
have used the latitude and longitude coordinates of each 
epicenter in relation to the origin of the coordinates, i.e., 
where both latitude and longitude are equal to zero (we 
have chosen the referential at the origin for simplicity). 
So, if a seismic event occurs with epicenter E with loca- 
tion (Oe-, 4>e), where 6e and (J)e are the values of latitude 
and longitude in radians of the epicenter, we are able 
to find the distances north-south and east-west between 
this point and the origin. These distances can be cal- 
culated, considering the spherical approximation for the 
Earth, by: 
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where 5"^* and 5'|f are, the north-south and east-west 
distances for the earthquake E, respectively, and R is the 
Earth radius, considered equal to 6.371 x 10^ km. With 
this computation we can identify the cell in the lattice 
for each event using the values of Se" and S]^ . Fig. fl] 
shows the globe divided into a square lattice where each 
square in the lattice is a potential node in the network. 

The seismic data used to build our network was 
taken from the Global Earthquake Catalog, provided by 
Advanced National Seismic System[T7], and it records 
events from the entire globe. The data spans all seismic 
events between the period from January 1, 1972 to De- 
cember 31, 2011. This catalog has a limitation because 
it is not consistent in all regions of the world; it includes 
events of all magnitudes for the United States of America 
but only events with m > 4.5 (in the Richter scale) for 
the rest of the world (unless they received specific infor- 
mation that the event was felt or caused damage) . There- 
fore, in order to obtain a more homogeneous distribution 
of data through the world, we have analyzed only events 
with m > 4.5 and we also excluded data that represent 
artificial seismic events ("quarry blasts"). In the end, we 
were left with 185 747 total events, where 82% of them 
happen near the surface of the world (depth < 100km). 



RESULTS 

Given the network build as described in Section [1 we 
have performed a few experiments to understand this 
structure. Following the procedure depicted in Fig. [T] 
the 185 747 yielded three different networks depending 
on the size used for the cells. Table |T] presents the sizes 
of the three networks. 



TABLE I. Three networks have been created from 185 747 
seismic events in our dataset. N represents the number of 
nodes and M represents the number of edges. Since we have 
a network constructed from consecutive seismic events, the 
number of edges, M, is always 1 less than the number of 
events. 



Network 



N 



M 



20 km X 20 km 65 355 185 746 
10 km x 10 km 104 516 185 746 
5 km x 5 km 144 974 185 746 



Scale- Free Property of the Seismic Network 

It has been shown the past by the work of Abe [71 [TS] 
that seismic networks for specific regions of the globe 
(e.g. California) appear to have scale-free properties, or 
in other words that the construction of the network em- 
ploys preferential attachment as described by [1] insofar 
that a node added to the network has a higher probabil- 
ity to be connected to an existing node that already has a 
large number of connections. This is somewhat trivial to 
understand because active sites in the world will tend to 
appear in the temporal sequence of seismic events many 
times. The preferential attachment states that the prob- 
ability P that a new node i will be linked to an existing 
node j, depends on the degree deg(j) of the node j, that 
is, P{i — ?► j) = deg(j)/^„deg(w). This rule generates 
a scale-free behavior whose connectivity distribution fol- 
lows a power-law with a negative exponent as shown in 
Eq.g 

In [71 [HI , earthquake networks were built for some spe- 
cific regions (California, Chile and Japan), and their con- 
nectivity distribution was defined as following a power- 
law. However, if we look carefully to the connectivity 
distribution and plot the cumulative probability of it, in- 
stead the probability density, we can observe that the 
power-law distributions that emerge from these network 
are truncated. According with [TH], there are two classes 
of factors that may affect the preferential attachment and 
consequently the scale-free degree distribution: the aging 
of the nodes and the cost of adding links to the nodes (or 
the limited capacity of a node). The aging effect means 
that even a highly connected node may, eventually, stop 



receiving new links as it normally occurs in scientific col- 
laboration networks [19] where scientists with time stop 
forming new collaborations maybe due to retirement or 
because they are already satisfied with the number of 
collaborators they have. The second factor that affects 
the preferential attachment occurs when the number of 
possible links attaching to a node is limited by physical 
factors or when this node has, for any reason, a limited 
capacity to receive connections, like in a network of world 
airports. Thus, these two factors impose a constraint to 
the preferential attachment and its power-law distribu- 
tion. When these factors are present, the distribution is 
better represented with a power law with an exponential 
cutoff, 



P(fc) 






(12) 



where a and fee are constants. 

In Fig. [2] we plot the cumulative probability distribu- 
tion for the earthquake network built for the Southern 
California (32°N - 37°N and 114°W - 122°W), using 
the data catalog provided by Advanced National Seismic 
System, where we considered all seisms with magnitude 
m > for the period between January 1, 2002 and De- 
cember 31, 2011. The total number of events are 147435. 
It is possible to observe in this plot that the data is bet- 
ter fit a power-law with exponential cutoff than a pure 
power law which is a good fit only for small values of 
k with an exponent 7 — 1 = 0.51, which is consistent 
with the value 7 = 1.5 reported in [7] for the probability 
density function. These results apply for a network built 
using cell sizes of 5 km x 5 km. It is noteworthy that in a 
probability density plot, the cutoff does not seem to exist, 
because the fluctuations are higher than in a cumulative 
probability plot. 



we verified if the magnitude distribution of seismic events 
in our data has the expected behavior. The Gutenberg- 
Richter law gives the rate of occurrence of earthquakes 
with magnitude larger or equal than m, 



F>{m) = 10" 
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(13) 



where, F> is the number of earthquakes with magnitude 
larger or equal than m and a and b are constants. 

As described in |S], this equation presents problems 
only on small values of magnitude. Since for the globe we 
are using data with m > 4.5, it is expected that our mag- 
nitude data have a good agreement with the Gutenberg- 
Richter law. This agreement is shown in Fig. [31 which 
gives the cumulative distribution of magnitudes. 
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FIG. 3. Log-linear plot for cumulative probability distribu- 
tion of the magnitude of earthquakes for the data used in 
this paper. The parameter b from the Gutenberg-Richter law 
(GR) has the value —1.08 ± 0.00058 and the associated value 
of the correlation coefficient, r, between the data (black filled 
circles) and the GR (red straight line) is 0.99903. 
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FIG. 2. Cumulative probability distribution of connectivity 
for the earthquake network in California using cell size 5 km x 
5 km. The solid lines represent two possible fittings: a power- 
law (black) and a power-law with exponential cutoff (red). 
The best fitting is for a power-law with exponential cutoff 
with a = -0.48 ± 0.0012 and kc = -279.0 ± 2.48. There are 
4 187 nodes in this network. 



Before constructing the network for the entire world. 



Looking at the world earthquake network constructed 
using the data described in Section [] we note that the 
aging-cost effect are visibly stronger in the connectivity 
distribution; the exponential cutoff is clearly visible in 
both the degree distribution and the cumulative degree 
distribution in Fig. |4J 

Fig. 4(a) represents the connectivity distribution for 
the global networks using the three different cell sizes for 
the global lattice. It is interesting to note that, com- 
paring these plots, we observe that the behavior is the 
same in all three cases, which indicates that the cell size 
does not change the complex features behind the global 
seismic phenomena. 

In Fig. |4(b)[ we have the same plot of Fig. 4(a) but 
using cumulative probability and only cell size 20 km x 
20 km. Note that the cumulative probability plot for the 
global network shows the same exponential cutoff behav- 
ior than for local network, as shown in Fig. [2J 

It is noteworthy that we did the same analysis in the 
connectivity distributions for the globe and for the Cal- 
ifornia using different thresholds for magnitude. The 
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FIG. 4. Connectivity distributions in the global earthquake 
network, (a) Plot for the cell sizes 20 km x 20 km (solid cir- 
cles), 10km X 10km (squares) and 5km x 5 km (cross), where 
the solid lines represent the best fit using power-law with ex- 
ponential cutoff. |(b)| Cumulative probability for the cell size 
20 km X 20 km. The solid lines represent a standard power- 
law and (black) a power-law with exponential cutoff (red). 
The best fit is the power-law with exponential cutoff with 
a = -0.96 ± 0.0011 and k^ = -17.0 ± 0.06. 



range of thresholds for magnitude were from 4.5 to 5.5 for 
the globe and from to 3.5 for Cahfornia. As expected, 
in all cases the distributions have the same behavior as 
shown in Fig. |4] 



To study these properties we need to introduce slight 
changes to our original network. The first is that the 
loops have to be removed, since we are looking for cor- 
relations between nodes and it only makes sense when 
these nodes are different. The second change is, that we 
move from a network with multi-graph characteristics to 
a weighted network. That is, if two nodes are linked by 
w edges in the original network, they will be linked by 
a single edge with weight w in the new version of the 
network. 

We have analyzed the seismic network for the entire 
world under two viewpoints: directed and undirected, 
where the cell size used in this construction was 20 km x 
20 km. The data used were the same as described in 
Section D Table Hi] shows the results obtained for the 
clustering coefficient (C) [20 and the average path length 



TABLE II. Results for the clustering coefficient (C) compared 
to the clustering coefficient of a random network of the same 
size (Crand) and the average path length {£) compared to the 
In A*', where A*' is 65 355 in network with cell size 20 km x 
20 km. 



Network 


C 


Crand 


e 


InN 


Directed 
Undirected 


0.007 
0.012 


4.2 x 10-^ 
4.2 X 10-^ 


17.19 
12.24 


11.08 
11.08 



From Table |TTj we note that both versions of the earth- 
quake network has small- world properties; the clustering 
coefficient is much higher than an equivalent for a random 
network, and the average path length has the same order 
of magnitude as the logarithm of the number of nodes. It 
is worth noticing that the regional earthquake networks 
built for California, Japan and Chile also are small-world 
[71 [TB] although the significance of small world at the 
global level is higher because we demonstrate long-range 
relations. 



Time interval between successive seismic events 



Small-world Property of the Global Seismic Network 

Small-world networks |10| have the general character- 
istic that they contain groups of near-cliques (dense ar- 
eas of connectivity) but long jumps between these areas. 
These two properties lead to a network in which the av- 
erage shortest path is very small and the clustering coef- 
ficient very high. 

Here we would like to test if the global seismic net- 
work has small-world properties. The consequence of 
such finding would be an indicative that seismic events 
around the world are correlated and not independent. 



We also studied the relationship between different re- 
gions of the world under another viewpoint, which is 
based on the analysis of the time intervals between the 
successive earthquakes in our global network. 

Previous studies have found that the probability distri- 
bution of time intervals between successive seismic events 
in small areas of the world (e.g. California and Japan) 
can be well described by nonextensive statistical mechan- 
ics [21 13] ■ We will verify in this section if these features 
are still present when we look to the entire world. 

In [3] , the authors use concepts from nonextensive sta- 
tistical mechanics to show that the cumulative probabil- 
ity distribution, P>, for the time interval between sue- 



cessive earthquakes in California and Japan follows a q- 
exponential, 



P>{At)^egi-pAt), 



(14) 



where, At is the time interval between successive events 
and /3 is a positive value. 

From Eqs. [9| and [T4| , it is possible to see that, if q > 1, 
when At ^ [/3(1 — q)]~^, the cumulative distribution 
represented in Eq. [M] approaches a power-law given by, 

p>(Ai)- Aii/(i-«r"^ 

Furthermore, from Eq. 
sides in Eq. [Ml 
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if we take the hiq in both 



lnq{P>{At)) = -13 At, 



(15) 



we can observe that the g-logarithmic of P>(Af) is linear 
with Ai with a slope — /3. 

Taking the worldwide data from Section [l we plot- 
ted the cumulative probability distribution for the time 
interval between successive earthquakes and we noticed 
that this distribution is well fitted by a g-exponential, 
indicating that the nonextensive behavior is also present 
when we look at seismic events from a global perspective. 
Fig. |5(a) shows the cumulative distribution on a log-log 
plot, where the histogram was made by using bins of 
size equal to 10 seconds. In Fig. |5(b)| we have the same 
distribution in a g log- linear plot, where the best value 
of q was found by analyzing the values of the correlation 
coefficients, as shown in the inset of Fig. [5"(b)| These cor- 
relation coefficients were taken in time intervals smaller 
than 51 000 seconds. 



CONCLUSIONS 

The use of networks to model and study relationships 
between seismic events has been used in the past for small 
areas of the globe. Here we demonstrate that similar 
techniques could also be used at the global level. More 
importantly, many of the techniques used in complex net- 
work analysis were used here to show that there seem to 
exist long-distance relations between seismic events. 

First we argued in favor of the long-distance relation 
hypothesis by showing that the network has small-world 
characteristics. Given the small-world characteristics of 
high clustering and low average path length, we were able 
to argue that seisms around the world appear not to be 
independent of each other. To strengthen this argument, 
we decided to do a temporal analysis of our network. 
Plotting the probability distribution for the time inter- 
vals between successive earthquakes, we have found that 
this distribution is well fitted by a g-exponential, indi- 
cating a behavior described by the non-extensive statis- 
tical mechanics, which obtain g-exponential distributions 
from the generalized Tsallis entropy. This non-extensive 
behavior also contributes to the long-distance relation 
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FIG. 5. Cumulative distribution of the time intervals between 
successive earthquakes with m > 4.5 in the entire world, (a 
log-log plot 



(b) qlog-lin plot, where the black dots represent 
the data and the red straight line represents the best fitting 
using the Eq. |15[ The slope of this plot gives /? = 1.43 x 
10~* ± 1.89 X 10~*s~^. Inset: linear correlation coefficient 
for each value of q, where we can see that the best fitting is 
reached by q = 1.08. 



hypothesis, since the non-extensive statistical has been 
used to explain many complex systems with long-range 
interactions and long-range temporal memory. Further- 
more, our results contribute for the conjecture of the con- 
nections between scale-free networks and non-extensive 
statistical mechanics, as proposed in [531 [53]. 

Another interesting approach we intend to do in the 
future relates to using community analysis or commu- 
nity detection to understand how seismic locations are 
grouped. We believe that given the long-range relations 
that we found here it is unlikely that the globe would be 
well organized around local communities of nodes (geo- 
graphical locations). 
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